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Abstract 

We confront the finite volume and small quark mass behaviour of the 
scalar condensate, determined numerically in quenched lattice QCD us- 
ing Neuberger fermions, with predictions of quenched chiral perturbation 
theory. We find that quenched chiral perturbation theory describes the 
numerical data well, allowing us to extract the infinite volume, chiral limit 
scalar condensate, up to a multiplicative renormalization constant. 
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Introduction 



Chiral symmetry breaking plays a central role in our comprehension of low energy 
QCD and understanding it from first principle calculations is of great importance. 
One of the cleanest ways of determining a condensate associated with the breaking 
of a global symmetry is through a finite-size scaling analysis. This technique 
has proved very successful in the study of scalar 0(N) models For chiral 
symmetry breaking in QCD, this would correspond to placing the system in a 
box and studying the scaling of the scalar condensate as a function of the volume 
V and of the quark mass m as the limit of restoration of chiral symmetry is 
approached {m ^ 0, V finite). 

The very small quark mass limit of QCD is expected to be well described by 
the lowest orders of chiral perturbation theory (xPT), which predict how the 
restoration of chiral symmetry takes place in a finite volume, as a function of the 
quark mass 0. The only free parameter entering the leading order contribution in 
the chiral expansion is the infinite volume quark condensate. Thus, a comparison 
of the mass and volume dependence of the finite volume quark condensate with the 
predictions of xPT provides a very powerful test of the hypothesis of spontaneous 
chiral symmetry breaking and permits an extraction of the infinite volume scalar 
condensate — S. Such a study requires, however, a good control over the chiral 
properties of the theory, which is difficult to achieve with traditional formulations 
of fermions on the lattice. 

The situation is different, however, when Dirac operators that satisfy the Ginsparg- 
Wilson (GW) relation Q are considered. Actions constructed from such opera- 
tors have been shown to have an exact lattice chiral symmetry [|^ . This symmetry 
ensures that the relations implied by chiral symmetry in the continuum, hold also 
on the lattice at finite lattice spacing a 0-0. For a review of the GW relation 
and its implications, we refer to |T^ . 



A particular realization of an operator satisfying the GW relation, which we will 
be using here, has been proposed by Neuberger 0: 

D^^[m + {l + s)il-j,Qm-'^% (1) 

where Q = Co75(l + s — -Dw)? -Dw is the Wilson Dirac operator, and the factor 
Co is a convenient normalization to keep the spectrum of bounded by 1. The 
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parameter s satisfies |s| < 1 and m is the bare quark mass. 

satisfies the GW relation at zero quark mass. In contrast to the standard 
Wilson formulation, the breaking of the chiral symmetry is soft, i.e. only due 
to the quark mass term. This opens the possibility to confront finite volume 
simulations with finite-size scaling predictions in the regime of restoration of the 
chiral symmetry. 

The complexity of the operator renders its numerical treatment very de- 
manding. We therefore restrict to the quenched approximation. The predictions 
of xPT must then be modified to take into account the effect of quenching. The 
finite size scaling of the quark condensate has recently been worked out using 
the framework of quenched chiral perturbation theory (q^PT) [jTl|. In particu- 
lar, analytical expressions for this scaling have been obtained in sectors of fixed 
topology. Operators satisfying the GW relation also satisfy an index theorem ||^. 
Thus, by computing the eigenvalues of at zero quark mass and identifying the 
zero modes, a clean separation of different topological sectors can be achieved, 
which is not possible with other formulations of lattice fermions. As we will see 
below, using the qxPT results in fixed topological sectors to interpret our numer- 
ical data proves very useful. A preliminary account of this work was presented 
at Lattice 99. 

Light quarks on a torus 

To study the volume dependence of the scalar condensate, we work on a four- 
dimensional torus of volume L^. Under the assumption that chiral symmetry 
is spontaneously broken, a description of QCD in terms of a chiral Lagrangian 
should be a good approximation at momenta p ^ AttF.^. To lowest order in 
p/AttFt^ and in the quark mass, this Lagrangian is given by 

£ = ^ Tr[ df,U^{x) d^U{x)] - S Re Tr[Me-'^/^^[/] (2) 

where U{x) = exp[i2n(x)/F^] G SU{Nf), Il{x) being the pion fields; M is the 
quark mass matrix, which we take to be proportional to the identity matrix (i.e. 
M = ml), and — S is the infinite volume and zero quark mass scalar condensate. 
In eq. (Q) we have included the expected 9 angle dependence. 
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Let us now consider the regime 



M^<1/L<F^, (3) 

where = 2mS/F^ to leading order in xPT. In this regime, the partition 
function is dominated by the zero mode of the U{x) field [0], since the action 
of the non-zero modes has a kinetic contribution that goes like F^L"^ ^ 1. The 
partition function then reduces, to leading order, to an integral over the SU (Nf) 
group manifold: 

Z = [ rf;7^eVEReTr[Me-«/-/t/o] ^ (4) 

JsU{Nf) 

where Uq is the global mode. We can also define the partition function restricted 
to fixed topology by Fourier transforming in 6* 0: 



Z, = [ ^ [ dUo e-'^" exp[V S ReTT[Me'^Uo]] 
Jo 27r Jsu{Nf) 

= [ dUo detiUoY exp[V S ReTr[Mf/o]]. (5) 

Ju{Nf) 

These integrals and their derivatives with respect to the quark mass have been 
known for a long time. For details see p|. 

In our case, however, we are interested in the quenched approximation. Recently 
a similar reasoning has been applied to quenched QCD. The main difference in 
the quenched case is that the chiral symmetry group is no longer SU{Nf)L x 
SU{Nf)ji X f/(l), but a graded Lie group U{1\1)l x t/(l|l)/j/f/A(l). According 
to the partition function for fixed topology is then given by 

Z^= [ c/[/o sdet([/o)'' exp[1/ S ResTr[M[/o]]. (6) 

Juaii) 



This integral has been computed analytically in terms of Bessel functions ||TT1 . 



By differentiating its logarithm with respect to the quark mass, m, the quark 
condensate for fixed topology is found to be 

= S z Mz)K{z) + I,+,{z)K_,{z)] + S^, (7) 

where z = m'EV and Ii^{z) , Ki,{z) are the modified Bessel functions. 
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This formula summarizes the scahng of the quark condensate in a periodic box 
with the volume and quark mass in the small mHV limit, as a function of only 
one non-perturbative parameter: the infinite volume condensate — S. For fixed 
volume, the limit as m — is given by 

S^=o = mJ:^V{l/2--f + log2-logmJ:V + 0{mJ:V\ogmJ:vf) 

S^=±i = + -mJ:^V (l + 0(mJ:V\ogmJ:V?) . (8) 

mV 2 ^ ^ 

where 7 is the Euler constant. These results have two interesting features that we 

wish to emphasize. First, there is a divergence ~ 1/m in sectors with topology. 

From the point of view of the underlying theory, this is not surprising since it 

corresponds to the contribution of the fermionic zero modes. Note however that 

these terms do not contain information about the infinite volume condensate and 

vanish in the infinite volume limit, as expected. The second interesting feature is 

the appearance of a logarithmic enhancement in T,,y=o, which is also peculiar to 

the quenched approximation. This term contains information about the infinite 

volume condensate. 

In principle, by fitting the dependence of the finite volume condensate in quark 
mass and volume to Monte Carlo data, we can extract the infinite volume conden- 
sate. However, the naive bare quark condensate that is measured on the lattice 
is UV-divergent. A simple dimensional analysis of the possible divergences shows 
that the bare scalar condensate has a leading cubic divergence. One important 
advantage of Neuberger's operator is that the coefficient of this leading divergence 
is known analytically. It is 6/(1 + s), for SU{3). The cubic divergence can then 
be subtrated exactly. However, after this trivial subtraction, the condensate is 
still divergent and has the form: 

ErV) - - (* *>. - (t^)t = C.'-^ + C:^ + E„. (9) 

1 + s a'^ a"^ a 

where m{a) is the bare lattice mass. The constants Ci are not known a priori 
and have to be determined, preferably non-perturbatively. The linear divergence 
proportional to m{af' is negligibly small for the values of the mass and the cutoff 
we consider in this work. However, the quadratic divergence is not and turns 
out to be very important numerically. The condensate extracted through a fit 
of the lattice data to eqs. and (|^), of course, still requires a multiplicative 
renormalization to eliminate a residual logarithmic UV divergence in Ej^. 
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After subtracting the unphysical contribution of fermionic zero modes to E,y, 
the finite volume condensate S^"^ vanishes, as expected, in the hmit of zero 



quark mass P, |T8[. Not surprisingly, the power divergences can be separated, 
in principle, from the physical contribution to the condensate, by a study of the 
volume dependence of S*"'', while keeping the quark mass small enough to stay 
in the region of validity of xPT. 

Even with the cubic divergence already subtracted, separating the physical con- 
densate from the remaining power divergences may not be easy, in practice, be- 
cause the statistical errors in these divergences can hide the small physical con- 
tribution. Thus, S*"* must be computed with very good accuracy. Clearly, the 
logarithmic enhacement of T,^=o in eq. could be very helpful in this respect; 
however, as we will see, extracting the condensate from the logarithmic term at 
zero topology requires much larger statistics than available to us at this time. We 
will concentrate instead on the study of the condensate in the topological charge 
one for minus one) sector. 



Numerical results 

For our numerical simulations we work on hypercubic lattices of size with 
periodic boundary conditions for both the gauge and the fermion fields. We work 
in the quenched approximation and use standard methods to obtain decorrelated 
gauge field configurations. 

In selecting the value of f3 = 6/ Qq, some care has to be taken. On the one hand, 
the quadratic divergence oc should not hide the physical effect. On the 

other hand, in choosing too small values of (3 there is the risk that Neuberger's 



operator falls into a different universality class [T^. Indeed, by computing the 
low-lying eigenvalues of Neuberger's operator at (3 = 5.7 and s = 0, we only 
found eigenvalues 0(1) and hence no light physical modes. A scan of the lowest 
eigenvalue of as a function of s showed that XmmiQ'^) decreased with increasing 
s, contrary to what is expected (and found) at larger values of f3. 
The situation at /? = 5.85 appeared to be different, however. The values of 
Amm(Q^) reach a maximum around s = 0.6, where the localization properties 
should also be optimal [|T2|. In accordance, the eigenvalues of Neuberger's opera- 
tor became very small, so that there is little doubt that light, physical modes are 
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present. Since the lattice spacing at /5 = 5.85 is a^^ ~ 1.5GeV |jl6[ we estimated 
that the quadratic divergence term would not hide the physical signal, at least 
for reasonable values of the physical condensate. 

A technical challenge is the numerical treatment of the square root appearing 
in Neuberger's operator. We have chosen a Chebyshev approximation for this 
task, which allows us to reach a well controlled accuracy. In order to avoid any 
systematic effects in the values of physical observables, we demand that 

\\X - Q^Pn,eiQyxf/i2\\X\\)' < 10-'' . (10) 



In eq. (|T0|) X denotes a random vector and P„^e denotes a standard Chebyshev 
approximation of the function in the range e < x < 1. Pn,e is a matrix- 

valued polynomial of degree n, which is constructed through numerically stable 
recursion relations |jl3|. We require tantamount accuracies for all inversions. We 
note in passing that with the requirement of eq. (0) also the GW relation itself 
is satisfied to a similar accuracy for zero mass. 

In order to decrease the degree of the polynomial employed, we have computed 
the 11 lowest eigenvalues of and their corresponding eigenvectors and have set 
e to be the value of the largest. The contributions of these lowest lying eigen- 
vectors are then treated exactly and projected out of the operator Q^. Through 
this procedure, near-zero modes of are taken into account automatically. All 
eigenvalue computations performed in our work are based on minimizing the Ritz 
functional ■ 



As pointed out in [|1^, it is advantageous for the computation of the eigenvalues 
of Neuberger's operator, and for its inversion as well, to stay in a given chiral 
subspace. This is possible because D]^Dj\f commutes with 75. 
We computed the scalar condensate at several values of the quark mass using a 
multiple mass solver on lattices of size 8^, 10^ and 12^. We checked through 
the calculation of the two lowest eigenvalues of Neuberger's operator to which 
topological sector each gauge field configuration belonged. We then obtained 
E™'' by computing 

^sub 1 — - I 1 1 



where the trace was performed in the chiral sector opposite to that with the zero 
modes and the gauge average was done in a sector of fixed topology z/. With this 
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definition, we take into account the contribution of all tlie non-zero eigenvalues 
of Dm to the condensate [} In this way, the term ~ 1/m in eq. (@) is absent. 
Three gaussian sources and standard inverters were used to compute the trace in 
eq. ([TT|) . Topological charge zero configurations are very rare at larger volumes. 
For this reason we did not compute the condensate in this sector on the 10^ and 
12^ lattices since the statistics we have gathered are too small. 
We show in fig. 1 our results for a^S^!f'^^/am on our lattice volumes as a function 
of bare quark mass. We have 15, 10 and 7 gauge configurations on our 8^, 10^ 
and 12^ lattices, respectively. The solid lines are a fit of the data for all volumes 
and masses to eqs. and (|^). This fit has only two parameters, namely the 
infinite volume, zero quark mass, scalar condensate — S and the coefficient of the 
quadratic divergence. We find a^S = 0.0032(4) and C2 = -0.914(8). 
Clearly, the formulae derived in xPT give a very good description of the numeri- 
cal data. The infinite volume condensate that we extract from this fit in physical 
units is — S(yU ~ 1.5 GeV) = — (221^g MeV)^, up to a multiplicative renormal- 
ization constant, which has not been computed yet for Neuberger's operator. We 
stress that the quoted error on the condensate is purely statistical. It does not 
include, for instance, the expected systematic errors from finite lattice spacing 
effects, nor the possible contributions from higher orders in chiral perturbation 
theory of 0{Ft^L)~'^ . An additional cautoniary remark is that the statistics for 
the largest volume, 12^, is rather small as indicated by the large statistical error. 
We plan to increase the statistics in the future and include in the analysis also 
higher topologies, which are more frequent at larger volumes. 
The condensate we obtain is quite close to that reported in ref. [|r^ using Wil- 
son fermions and a different method. However, a meaningful comparison can 
only be made when our systematic errors are quantified and the multiplicative 
renormalization included. 

Recently the authors of |TB[ also studied the quark condensate as a function of 
quark mass and volume, using Neuberger's operator. However, a comparison with 
the predictions of qxPT in fixed topological sectors was not attempted and no 
definite conclusion on the existence or value of the infinite volume condensate 
was reached. 



^With this definition the real eigenvalues at the cut-off level, m + 2/a, are doubly counted. 
Although it is a completely negligible effect, we took it into account. 
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Figure 1: We show the quark mass dependence of the scalar condensate in three 
volumes, 8"^ (circles), 10^ (squares) and 12^ (triangles). The solid curves represent 
a fit of the data to eqs. (|^) and (|^). 
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According to Random Matrix Theory (RMT), the value of S may also be ex- 
tracted from the distribution of the lowest non-zero eigenvalue Amin(-DN), defined 
by the square root of the lowest non-zero eigenvalue of D^jD^ at zero quark mass. 
We have only gathered a reasonable statistics for the smaller lattice of size 8^. In 
the topological charge z/ = 0, 1 sectors, the corresponding distributions are given 
by p||: 



P.=o{z) = \e'-^'\ Pu=±i{z) = '-h{z) e-3^^ (12) 

where z = X^i^{D]<()T,V . Recently, the authors of ref. found very good 
agreement with these distributions in very small lattices. Inserting our value 
of a^S = 0.0032(4) in the distribution for zero topology of eq. (|T2p, we get for 
the expectation value of this eigenvalue (Amin(-DN)) = 0.135(15) (where the error 
comes from the statistical error in the condensate), while from our data on the 8^ 
lattice (with 41 topology zero configurations) we obtain {Xmm{DN)) = 0.170(12). 
In topology one sectors, the expected value for the (Amin(-DN)) = 0.237(17). From 
the data, again obtained on the 8^ lattice ,with an accummulated statistics of 29 
configurations, we obtain {\min{DN)) = 0.218(13). In addition, we generated 
a sample of eigenvalues according to the distribution of eq. (|T2|) for z/ = ±1 
with a^S = 0.0032(4) and a statistics identical to that of the corresponding 
simulation. The resulting mean value of Amin(-DN) and error are fully compatible 
with those given by our simulation, indicating that our data do not suffer from 
autocorrelation effects. This provides a nice cross check on the value of S obtained 
from our finite-size scaling analysis. 

We finally briefly comment on the difficulty in measuring the condensate in the 
topology zero sector. This is due to the logarithmic enhancement, which can be 
shown to originate from the contribution of the single lowest eigenvalue, Ainin(-DAr) 
to the condensate, if the distribution of this eigenvalue is that given by RMT in 
eq. ([T^). As is clear from eq. (|T2|), the distribution of the lowest eigenvalue 
does not have a gap and it is easy to check that the contribution of this single 
eigenvalue to the condensate, Smin, has a logarithmic IR divergence in the sector 
of zero topology: 

1 r 2m 

Smin = T7 dz P,=o{z) , , , = -mSVlogmSl^ + O(mSV), (13) 
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which reproduces exactly the logarithmic dependence in eq. (|^). We have gener- 
ated a sample of eigenvalues according to eq. (|T^) with a^S = 0.0032. In doing 
so we find that reconstructing the logarithmic behaviour of the u = condensate 
requires a statistics much larger than the one available to us. This is refiected in 
the results from our actual simulation where we find that the = condensate 
on the S'* lattice displays very large statistical errors. For this reason, we have 
not included these data in our determination of the scalar condensate. 
Conclusion 

Chiral perturbation theory assumes that chiral symmetry is spontaneously bro- 
ken in QCD. Under this assumption it provides, for small enough values of the 
quark mass and large enough volumes, the mass and volume behaviour of the 
scalar condensate. This behaviour is determined to lowest order by only one 
free parameter, namely — S, the scalar condensate in infinite volume and for zero 
quark mass. 

Recent developments in lattice QCD have revealed that, contrary to a long- 
standing belief, chiral symmetry can be realized on the lattice. This theoretical 
advance is connected to the Ginsparg-Wilson relation. Neuberger proposed a 
particular operator that satisfies this relation and we have used this operator in 
our numerical work. 

Although Neuberger's operator is very difficult to treat numerically, it can be 
used in practice. In this work we computed the scalar condensate on lattices of 
various sizes and for a number of quark masses, in the regime of chiral symme- 
try restoration. The results of this numerical computation are shown in fig. 1, 
where we confront our numerical data with the finite volume and mass behaviour 
predicted by quenched chiral perturbation theory. Obviously, chiral perturbation 
theory describes the numerical data well, providing evidence for the spontaneous 
breaking of chiral symmetry. 

Although our results are very encouraging, a number of cautionary remarks have 
to be made. The lattices we used are rather small and it would be desirable to 
probe the system further on larger lattices. In addition, it would be important 
to repeat the calculation at a larger value of /3 to estimate the lattice spacing 
effects. Here we were only able to determine a value of the scalar condensate up 
to a multiplicative renormalization constant, which would clearly be needed for 
quoting a physical value. Finally, all our results are obtained in the quenched 
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approximation but, given the complexity of Neuberger's operator, it would be 
very difficult to go beyond this approximation. 



During the completion of this work, a paper using Neuberger's operator to 
compute the scalar condensate appeared. The data presented in this paper are, 
however, taken in the strong coupling regime on only one lattice (with small 
L/a=4) and can hence not directly be compared to our work. 
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